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The a 4 contribution to the lepton g — 2 from a gauge-invariant set of 18 Feynman diagrams 
containing a light-by-light scattering subdiagram internally has been reevaluated by a method inde- 
pendent of the previous approach. Comparison of two methods revealed a program error in the first 
version. Correcting this error, the contributions of these 18 diagrams become -0.990 72 (10)(a/7r) 4 
and -4.432 43 (58)(a/-7r) 4 for the electron and muon g — 2, respectively. The correction is not large 
enough to affect the comparison between theory and experiment for the muon g — 2, but it does 
alter the inferred value for the fine structure constant a -1 by 6 ppb. 

PACS numbers: 13.40.Em,14.60.Cd, 14.60.Ef, 12.20.Ds 



Precise theoretical and experimental values of lepton 
anomalous magnetic moments (a/ = (<?; — 2)/2) provide 
one of the most stringent tests of QED [Q . In case of the 
electron (positron) , the experimental value did reach the 
precision of 4.3 ppb Currently, the theoretical uncer- 
tainty is dominated by that of the fine structure constant 
a. The most precise a available now is from the atom in- 
terferometry experiment which has 7.4 ppb precision. 
Since the muon anomalous magnetic moment a M is sen- 
sitive to short-distance physics, high precision measure- 
ment (1.3 ppm) of a M at Brookhaven National Laboratory 
may be able to open the first window to "new physics" 
j|. Before taking the discovery of "new physics" in the 
muon g — 2 seriously, however, we must make sure that 
the old physics, namely the Standard Model, is known 
with sufficiently high precision. 

The largest source of theoretical uncertainty (0.7 ppm) 
for a ^ is the hadronic contribution j^, Unfortunately 
we are currently unable to deal with the hadronic cor- 
rection from first principles because of nonperturbative 
nature of QCD. On the other hand, the QED correction 
can be treated precisely by perturbation theory. In order 
to achieve the precision comparable to that of measure- 
ments, the QED calculation for lepton g — 2 must include 
terms of up to eighth-order of perturbation theory. Lead- 
ing contributions of tenth-order are also relevant for 

The purpose of this paper is to correct a program error 
in the previous calculation of gauge-invariant set of 18 
Feynman diagrams contributing to the a 4 QED term . 
This was accomplished by constructing alternative forms 
of integrals for these diagrams. As a consequence all 
891 Feynman diagrams contributing to the eighth-order 
term of a R , and additional diagrams contributing to a M — 



order of magnitude improvement in numerical precision 
of all a 4 terms of Ou and a e . The results will be reported 
shortly elsewhere M |To[ | . 

The contribution of the QED diagrams to a M can be 
written in the general form 

a M (QED) = Ax + A 2 {m tl /m e ) + A 2 {m ll /m T ) 

+ A 3 (m tl /m e ,m tJ ./m T ), (1) 

where m e , to m , and m T are the masses of the elec- 
tron, muon, and tau, respectively. A similar equa- 
tion holds for a e . Throughout this article we shall 
use the values m e = 0.510 998 902 (21) MeV/c 2 , m M 
=105.658 357 (5) MeV/c 2 , and m T = 1 776.99 (+29,- 
26) MeV/c 2 JuJ . 

The renormalizability of QED guarantees that the 
functions A\, A 2 , and A3 can be expanded in power series 
in a/ir with finite calculable coefficients: 
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.... i = 1,2,3. 

(2) 

Af\ a{ 4) and Af ] are known analytically O. Most 
terms contributing to A\ have not yet been obtained 
by analytic means. The current uncertainty in the value 
of A\ is a consequence of the fact that, at present, it 
must be obtained by numerical integration. Its precision 
is being improved by an extensive computer calculation 
right now ||. For the purpose of evaluating a AI (QED), 
however, it is sufficient to use A\ derived from the mea- 
sured value of the electron anomaly a e fl3ll corrected 
for small contributions of muon, hadron, and weak in- 



teractions. The terms A 2 ^ (m^jm^ and A^' '(m M /m e ) 
have now been verified by independent calculations ar e known exactly Q. The situation is quite different 

for A 2 (m,ij,/m e ) since most terms contributing to it are 



(6), 



and/or checked by analytic comparison with lower-order 
integrals. This enables us to pursue with confidence an 



known only by numerical means. 
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FIG. 1: Vertex diagrams containing a light-by-light scatter- 
ing subdiagram internally. There are altogether 18 such dia- 
grams. 



There are altogether 469 Feynman diagrams contribut- 

(8) 

ing to A 2 (m M /m e ). Of these diagrams 343 have been 
checked by more than one independent methods, some 
of which being even analytic. Other 108 diagrams, all of 
which contain an external light-by-light scattering sub- 
diagram, have been checked analytically by comparison 
with the exactly known sixth-order vertices |l^] . Unfor- 
tunately, this was not the case for the remaining 18, all 
generated by inserting a light-by-light scattering subdi- 
agram internally in a fourth-order vertex diagram (see 
Fig. 1). Here external and internal means whether one 
of the attached photon lines represents an external mag- 
netic field or not. 

Let us now describe how the error was discovered and 
corrected. These 18 diagrams form a gauge-invariant set 
and share the same basic algebraic structure. Unfortu- 
nately, it was not possible to examine their analytic struc- 
ture by comparison with lower-order diagrams since they 
are not reducible to such diagrams in the UV and/or IR 
limits. Besides having been checked by two people work- 
ing independently |p^| , the only check made was mutual 
consistency among these 18 diagrams. This is not suffi- 
cient to eliminate the possibility that they share the same 
program error. Clearly, in order to enhance the credibil- 
ity of the QED calculation of ai , it is highly desirable to 
reevaluate these 18 diagrams by more than one method. 

This paper reports the consequence of our effort to re- 
examine the previous result by construction of new and 
independent integrals. Before describing the new ap- 
proach, let us quickly go over the initial approach. In that 
approach (which will be called Version A) jl5| , we put to- 
gether vertex diagrams, for instance llj(l),llj(2),llj(3), 
all obtained from the self-energy-like diagram 11 j of Fig. 
2 as terms linear in the magnetic field in the weak field 
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FIG. 2: Self-energy-like diagrams in which lepton lines prop- 
agate in the magnetic field. First-order terms in the weak 
magnetic field expansion correspond to the diagrams of Fig. 





expansion. This is to take advantage of their shared 
structure and the tendency of partial cancellation among 
them. With the help of Ward-Takahashi identity the 
sum A" (p, q) of all vertex diagrams thus related to a self- 
energy diagram S(p) can be expressed as 



A"(p,«) 



9A M (p, q) 



q=0 



0S(p) 
dpu 



(3) 



where (p + q/2) 2 = (p — q/2) 2 = mf. The g — 2 term is 
projected out from the right-hand side of (|^) . In terms of 
Feynman parameters zi, z%, ■ ■ ■ Zn, the n-th order mag- 
netic moment derived from (0) has a form 



M (2n) 



(^)V-l)! J{dz){ 



E + C 1 



n - 1 UW 



2i/n-l 



(4) 



where bold letters E, C, N, and Z stand for parts of the 
projection operator adapted to the first and second terms 
on the right-hand-side of (||) . U is the Jacobian of trans- 
formation from momentum-space variables to Feynman 
parameters. V~ l is obtained by combining all propaga- 
tors into one with the help of Feynman parameters. See 
p6[ for precise definitions of projection operators, U, V 
and (dz). 

In order to check the validity of Version A based 
on Fig. 2 we evaluated these diagrams by a second 
method, called Version B, which is actually a straight- 
forward parametrization of individual vertex diagrams of 
Fig. 1, without relying on the Ward-Takahashi identity. 
In this approach it is convenient to evaluate the sum 
jkl{n) = llj(n) + llk(n) + lll(ri), where n = 1,2, 3, since 
partial cancellation of singular terms occurs resulting in 
less singular behavior. 

Diagrams of Fig. 1 (or Fig. 2) form a (formal) gauge- 
invariant set. But individual diagrams are UV-divergent, 
and must be regularized in advance, for instance by di- 
mensional regularization, to enforce gauge invariance. 
For numerical evaluation, however, it turns out to be 
more convenient to combine it with a subtractive reg- 
ularization. Let F m (d) be one of the integrals defined 
in (^-dimension, where m is any one of . . . , 111(3). 

Let G m (d) be the subtraction term containing the light- 
by-light scattering tensor with zero external momenta 
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n AlI y (T p(0, 0, 0, 0) as well as terms containing a sixth-order 
charge renormalization diagram. Let us rewrite F m (d) 
symbolically as 

(F m (d) - G m {d)) + G m (d), (5) 

where "symbolically" means that subtraction is per- 
formed on the integrand before the integration is carried 
out. Now, we can safely take the limit d — » 4 for the first 
term since its integrand does not cause UV divergence. 
Of course, the second term G m (d) is singular for d —> 4. 
However, gauge invariance guarantees that the sum of 
G m (d) over all diagrams of Fig. 1 vanishes for any value 
of dimension d: 

J2Gm(d)=0. (6) 

m 

Thus we have only to compute (F m (4) — G m (4)) in the 
end. 

The integrands of both Version A and Version B were 
generated by an algebraic program FORM |l7j . Numer- 
ical integration is carried out by an adaptive-iterative 
Monte-Carlo routine VEGAS g§. 

When we compared the numerical results of Version 
A and Version B, we were surprised to find that their 
values were significantly different. After an extensive de- 
tective work, we located a programming error in Version 
A, which resulted from an incomplete implementation of 
the E operation of (|4|) in the algebraic manipulation pro- 
gram: It left out some terms referring to the light- by- light 
loop subdiagram. Such a program error can be readily 
detected if the integral exhibits UV or IR divergence after 
renormalization is carried out. Unfortunately the partic- 
ular error in the 18 diagrams caused no divergence and 
escaped scrutiny of two people. Once this error was cor- 
rected, both approaches gave identical numerical results. 
Thus we now have two sets of independent codes for the 
18 diagrams that have been fully verified. 

Numerical evaluation of these diagrams requires an 
enormous amount of computational effort. A systematic 
algorithm of computation to minimize human error is in- 
dispensable. Such a scheme was developed originally for 
the calculation of the sixth-order lepton g — 2 [[ll| [HJ, and 
were thoroughly tested over the years |^. It was later 
extended to the eighth-order pi]] . Early results for the 
eighth-order term remained rather crude for many years. 
This is mainly due to the enormous size of the integrands 
which could not be handled adequately by the computers 
then available. More precise values have become avail- 
able only in this decade thanks to the development of 
massively parallel computer, which enabled us to vastly 
increase the sampling statistics of VEGAS. 

Enlarging sampling statistics, however, amplified the 
difficulty caused by a previously poorly understood prob- 
lem in estimating errors in a computer calculation. This 
arises from the fact that computer calculation always 



TABLE I: Muon g — 2 contributions from the diagrams of Fig. 
1. In Version A the Ward-Takahashi-summed llj = llj(l) + 
llj (2) + llj(3), Ilk = ZZfc(l) + llk{2) + llk(3), and III = 111(1) + 
111(2) +111(3) are calculated, while in Version B jkl(l,3) = 
llj(l) + llj(3) + ZZfe(l) + llk(3) + 111(1) + 111(3), and jkl(2) = 
ZZj'(2) + llk(2) + 111(2) are calculated. 





Version A 




Version B 


llj 


6.389 802 (460) 


jfcZ(l,3) 


-3.509 978 (802) 


Ilk 


-7.763 474 (537) 


jkl(2) 


-0.921 589 ( 89) 


III 


-3.059 704 (452) 






sum 


—4.433 376 (840) 


sum 


-4.431 567 (806) 



TABLE II: Electron g — 2 contributions from the diagrams 
of Fig. 1. In Version A the Ward-Takahashi-summed 
llj = llj(l) + llj(2)+llj(3), Ilk = llk(l) + llk(2) + llk(3), and 
III = 111(1) + 111(2) + 111(3) are calculated, while in Version B 
jkl(l,3) = llj(l)+llj(3) + llk(l) + llk(3) + 111(1) + 111(3), and 
jkl(2) ee llj (2) + llk(2) + 111(2) are calculated. 





Version A 




Version B 


llj 


2.551 223 (78) 


jkl(l,3) 


-0.872 717 (138) 


Ilk 


-1.873 801 (72) 


jkl(2) 


-0.117 959 ( 28) 


III 


-1.668 182 (80) 






sum 


-0.990 760 (133) 


sum 


-0.990 675 (141) 



deals with a finite number of digits. This means that 
the error estimate based on the assumption of normal 
distribution of errors must be modified to take the effect 
of rounding-off of digits into account. In our formulation 
in which subtractive renormalization of UV divergence as 
well as separation of IR divergences are carried out on the 
computer, this digit deficiency problem can distort error 
estimates seriously, or even prevent further iteration. Be- 
sides increasing the number of effective digits from real*8 
to real* 16 arithmetic, which was the most obvious and 
effective cure, various methods had to be devised to deal 
with the digit deficiency error (22). 

New results of numerical integration of by Versions 
A and B arc listed in Table I. The values of llj, Ilk, and 
III listed were obtained using 5 x 10 9 sampling points per 
iteration and iterated 110, 219, 220 times, respectively. 
They were evaluated on vl cluster at Cornell Theory Cen- 
ter. The calculation of Version B was carried out on 
Fujitsu VPP700E at the Computer and Information Di- 
vision of RIKEN. For jfc/(l,3) 4.6 x 10 9 sampling points 
per iteration were used for 131 iterations. The program 
jkl(2) shows less singular behavior. It was evaluated us- 
ing 4.6 x 10 9 sampling points per iteration and iterated 
60 times. 

Results for the electron are listed in Table II. The val- 
ues of llj, Ilk, and III listed were obtained using 2 x 10 9 
sampling points per iteration and iterated 160, 220, 180 
times, respectively, on vl cluster at Cornell Theory Cen- 
ter. For jkl(l,3) 4.6 x 10 9 sampling points per itera- 
tion were used for 60 iterations on VPP700E. For jkl(2) 
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4.6 x 10 9 sampling points per iteration were used and 
iterated 60 times. 

Combining the results from Tables I and II, which we 
treat as statistically independent, we obtain the best es- 
timate of the contribution to a M from the 18 Feynman 
diagrams of Fig. 1: 



n{ S) ,f„ = -4.432 43 (R8) 



"IV(d) 



jjx 



m 

[3] 

[4] 
[5] 



and a corresponding result for a e 



R. S. Van Dyck, Jr., P. B. Schwinberg, and H. G. 
Dehmelt, Phys. Rev. Lett. 59, 26 (1987). 
A. Wicht et al, in Proc. of 6th Symposium on Frequency 
Standards and Metrology edited by P. Gill (World Sci- 
entific, Singapore, 2002), pp. 193 - 212. 
G. W. Bennett, et al. (Muon (g-2) Collaboration), Phys. 
Rev. Lett. 89, 101804 (2002). 

For the hadronic vacuum polarization contribution , see 
M. Davier, S. Eidelman, A. Hocker, and Z. Zhang, hej> 
ph/0208177| ; K. Hagiwara, A. D Martin, D. Nomura, 
and T. Teubner, |hep-ph/0209187|; F. Jegerlehner, talk 



,(8)e 
l IV(d) 



-0.990 72 



(8) 



Here, superscripts (8)fi and (8)e refer to the eight-order 
muon and electron g — 2, the subscript IV to the group 
of all diagrams containing light-by-light scattering sub- 
Cams, and (<7) to its snhgronp shown in Fig 1 which 



at workshop " Hadronic Contributions to Anomalous 
Magnetic Moment of the Muon", Centre de Physique 
Theorique, Marseille, 14 -16, March 2002. 
[6] For the hadronic light-by-light contribution, see M. 
Knecht and A. Nyffeler, Phys. Rev. D 65, 073034 (2002); 
M. Knecht, A. Nyffeler, M. Perrottet and E. de Rafael, 
Phys. Rev. Lett. 88, 071802 (2002); M. Havakawa and 



diagx 



consists nf all diagrams containing an internal light-hy- 



T. Kinoshit a, Phys. Rev. D 66, 019902 (2002) and |hep- 
ph/0112102]; J. Bijnens, E. Pallante and J. Prades, Nucl 



light scattering diagram. 

The new results (R) and (||) supersede the earlier values 
-3.4387 (533) (q/tt) 4 § and -0.7503 (60) (07V) 4 H, 



Phys. B626, 410 (2002); I. Blokland, A. Czarnecki and 
K. Melnikov, Phys. Rev. Lett. 88, 071803 (2002); M. 
J. Ramsey-Musolf and M. B. Wise, Phys. Re v. Lett. 



respectively. The effect of this modification on a M (QED) 
is less than 1% of the overall eighth-order term (which 
is of the order of 130 (a/ir) 4 ), and thus does not af- 
fect comparison of experiment and theory significantly. 
On the other hand, the effect on a e is ~ —7.0 x 10~ 12 
which is about 16 % of the entire eighth-order term and 
is larger than the measurement uncertainty 4.3 x 10 -12 . 
As a consequence, it reduces the inverse fine structure 
constant oT x obtained from theory and measurement of 
a e by ~ 0.82 x 10~ 6 or ~ 6 ppb. Currently all a 4 terms 
are being upgraded by an extensive numerical integra- 
tion. Precise values of a M (QED) and a e (QED) including 
these terms will be reported in [l(J . 
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